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Abstract 

In experiment, the multiplicity distributions of inelastic processes 
are truncated due to finite energy, insufficient statistics or special 
choice of events. It is shown that the moments of such truncated 
multiplicity distributions possess some typical features. In particular, 
the oscillations of cumulant moments at high ranks and their negative 
values at the second rank can be considered as ones most indicative 
on specifics of these distributions. They allow to distinguish between 
distributions of different type. 

1 Introduction 

Studies of multiplicity distributions of high-energy inelastic processes have 
produced many important and sometimes unexpected results (for the reviews, 
see, e.g., |], ^, 3). The completely new region of very high multiplicities will 
be opened with the advent of RHIC, LHC and TESLA accelerators. 

Theoretical approaches to multiplicity distributions in high-energy pro- 
cesses have usually to deal with analytic expressions at (pre) asymptotic en- 
ergies which only approximately account for the energy- momentum conser- 
vation laws or with purely phenomenological expressions of the probability 
theory. The multiplicity range extends in this case from zero to infinity. 

In experiment, however, one has to consider distributions truncated at 
some multiplicity values in one or another way. These cuts could appear 
due to energy limitations, low statistics of experimental data or because of 
special conditions of an experiment. Energy limitations always impose the 
upper cotoff on the tail of the multiplicity distributions. Low statistics of data 
can truncate these distributions from both ends if it is insufficient to detect 
rare events with very low and/or very high multiplicity. Similar truncations 
appear in some specially designed experiments when events within some 
definite range of multiplicities have been chosen. 

It would be desirable even in these cases to compare the distributions 
within those limited regions with underlying theoretical distributions. The 
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straightforward fits are sometimes not accurate enough to distinguish be- 
tween various possibilities because the probability values vary by many orders 
of the magnitude. More rigorous approach is to compare different moments 
of the truncated distributions. Among them, the cumulant moments K q seem 
to be most sensitive to slight variations (and, especially, cuts and shoulders) 
of the distributions. They often reveal such tiny details of the distributions 
which otherwise are hard to notice. 

In particular, QCD predicts quite peculiar behaviour of cumulant mo- 
ments as functions of their rank q. According to solutions of the equations 
for the generating functions of the multiplicity distributions in the asymp- 
totic energy region, the ratio of cumulant moments K q to factorial moments 
F q usually denoted as H q = K q /F q behaves as q~ 2 and at preasymptotic 
energy values reveals the minimum || at q ~ 5 with subsequent oscillations 
at higher ranks |7J. Such a behaviour has been found in experiment at 
presently available energies || |9| . The solutions of the corresponding equa- 
tions for the fixed coupling QCD also indicate on similar oscillations [[IIJ. At 
asymptotics, the oscillations should disappear and H q becomes a smoothly 
decreasing and positively definite function of q, as mentioned above. 

Neither of the distributions of the probability theory possesses these fea- 
tures. Among them, the negative binomial distribution (NBD) happens to 
be one of the most successful ones in the description of global features of 



the multiplicity distributions [|Tl|] . Let us remind that the negative binomial 
distribution is defined as 

T{n + k) a n (l + ay n - k , (1) 



T(n + l)T(k) 



where a = (n)/k, (n) is the mean multiplicity, k is an adjustable parameter, 
and the normalization condition reads 



E^n = L (2) 

n=0 



Its generating function is 



G(z) = ±P n (l + zT=(l- Z -^-) . (3) 



n=0 
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The integer rank factorial and cumulant moments, and their ratio are 



1 dPG{z) T(q + k) 

q (n)* dz* U= ° V{k)ki ' 1 ] 

1 <PlnG(*) , T(g) 

r( g )r(fc + i) 
9 = r(fc + g ) • (6) 

ifg-moments at the parameter k = 2 behave as 2/g(g+l), i.e. with the power- 
law decrease reminding at large q that of QCD, however, with a different 
weight factor. Therefore, at first sight, it could be considered as a reasonably 
good analytic model for asymptotic behaviour of multiplicity distributions. 
It has been proclaimed fi"2 , 13, 14, 15] that the superposition of two NBDs 



with different parameters and their cutoff at high multiplicities can give rise 
to oscillations of H q and better fits of experimental data at preasymptotic 
energies. Nevertheless, the fits have not been perfect enough. 

Let us compare first the asymptotic QCD predictions with NBD fits at 
different values of the adjustable parameter k. The values of D q = q 2 H q are 
plotted in Fig. 1 as functions of q for the asymptotic QCD (where they are 
identically equal to 1) and for NBD at different k. At k = 2, they exceed 1 
tending to 2 at large q. At larger values of k, all D q are less than 1 except 
D% — 1 at k — 3. Surely, the identity D\ = 1 is valid for any k due to the 
normalization condition. 

To get asymptotic QCD results with all D q = 1 from the expressions 
similar to NBD, one would need to modify NBD in such a way that the 
parameter k becomes a function of n. Thus some effective values of k should 
be used to get QCD moments D q — 1 at various q. These effective values are 
shown in Fig. 1 at the corresponding QCD points. They are obtained as the 
solutions of the equation 



q-l 



k(q) 



n [ 1+ nj =q2 > ^ 



n=l 



which follows from Eq. (|6|) for H q = q 2 . They show that k somewhat 
decreases from 3 to some values exceeding 2 with increase of q. This reflects 
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the well known fact that the tails of distributions are underestimated in NBD- 
fits ]14| compared to experimental data in the preasymptotic region. Also, 
the amplitude of oscillations and their periodicity are not well reproduced by 



a single truncated NBD ||14]| , and one has to use the sum of at least two NBDs 
to get a better fit. However, rather large values of k were obtained in these 
fits. It implies, in fact, that the fit is done with the help of two distributions 
very close to Poissonian shapes because the Poisson distribution is obtained 
from NBD in the limit k — > oo. Therefore, the tails are suppressed very 
strongly. 

Here, we will not try to fit experimental data focusing our efforts on 
qualitative changes of moments when NBD is truncated, especially, as applied 
to studies of very high multiplicities. 

In QCD considerations based on the equations for the generating func- 
tions for quark and gluon jets, the preasymptotic (next-to leading order etc) 
corrections give rise to oscillations of H q . Even though they are of the higher 
order in the coupling strength, they appear mainly due to account of energy 
conservation in the vertices of Feynman diagrams but not due to considering 
the higher order diagrams which are summed in the modified perturbation 
theory series (see 0). In the phenomenological approach, this would effec- 
tively correspond to the cutoff of the multiplicity distribution at some large 
multiplicity. Therefore, we intend here to study how strongly such a cutoff 
influences the NBD-moments, whether it produces oscillations of the cumu- 
lant moments, how strong they are, and, as a more general case, consider 
the moments of NBD truncated both at low and high multiplicities. This 
would help answer the question if the shape of the distribution in the lim- 
ited region can be accurately restored from the behaviour of its moments. 
It could become especially helpful if only events with very high multiplic- 
ities are considered in a given experiment because of the above mentioned 
underestimation of tails in the NBD-fits. 



2 Truncated NBD and its moments 

In real situations, the multiplicity distribution is sometimes measured in some 
interval of multiplicities and one can try to fit by NBD the data available only 
in the restricted multiplicity range. Therefore, we shall consider the negative 
binomial distribution within the interval of multiplicities m < n < N called 
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P^ and normalized to 1 so that 

N 

E = i- (8) 

n=m 

Moreover, due to above reasoning and to simplify formulas we consider here 
only the case of k = 2. The generalization to arbitrary values of k is straight- 
forward. 

The generating function of the truncated distribution G c (z) can be easily 
found as 

N f(z) 
Gc{z) = E m Pi c) (l + z) n = G(z)(l + zT J j^ y (9) 

where 

/(*) = 1 + m{\ - x) - [1 + (N + 1)(1 - a;)]^-"^ 1 , (10) 
x = b{l + z), b = -^-. (11) 

J. — r~ (X 

Correspondingly, 

/(0) = /(s = 0) = /(* = &). (12) 

Using the above formulas for the factorial moments, one gets the following 
formula for the moments of the truncated distribution expressed in terms of 
the NBD-moments (HI): 



r(g + i) 

_ o r(r + l)r(g-r + l) 



(n) q F (c) = V V ; (n) q ~ r F 



rg)^ - b) + 1 ~ r] - r^Sy^M^ + 1)(1 - 6) + 1 - r] 

m(l - b) + 1 - & JV ~ m+1 [(TV + 1)(1 - 6) + 1] ' 1 j 

where (n) c is the mean multiplicity of the truncated distribution. It is related 
to the mean multiplicity (n) of the original distribution as 

(1 - b)[(N +1)(N + 2 )b N - m+1 - mjm + 1)] 
W Wc 1 + m{\ - b) + b N - m+1 [{N + 1)6 - N - 2] ' 1 J 



Inserting formula (§]) in fll3|), one gets 

/ w F (e) = (g-r + l)r(g+l) 
Wc 9 h 2^T(r + l) {U) 



r(m+2-r) L'"^ ~~ "V T - 1 - ~~ ' J ~~ r(7V+3-r) ^ 



^±^r m (i _ b ) + i _ r ] _ r(jv+3) ^-m+i^jy + i)(i _ 6) + i _ r ] 



m(l - 6) + 1 - 6 iV ~ ni+1 [(TV + 1)(1 - 6) + 1] 

5 



(15) 



For computing it is more convenient to use the formula fll5"|) in the following 
form: 




q a 



M 1 + 7(0)(?TT)5- ( « + 1 - r)x 



r=l 



[(— hi- r)8(m + 1 - r) JJ(m + 2-i) 



-(^I + 1 _ r ) & iV-m+l pj (]V + 3 _ •)] 1 (16) 
a + 1 i=l J 

These expressions can be used also for the distributions truncated at one 
side by setting m = or iV = oo. 

The cumulant moments can be calculated after the factorial moments are 
known from Eq. ([13]) according to the identities 

ql r (g) 
i w I TvFT \Kq_ m F m . (17) 

This formula is a simple relation between the derivatives of a function and 
of its logarithm (see Eqs @ and (|5])). Therefore it is valid for both original 
and truncated distributions. 

For the Poisson distribution, the ratios H q are identically equal to zero, 
and are given by Eq. (|[) for NBD while truncation induces new features. In 
Figures we show the behaviour of the ratios H q as functions of the rank q for 
the truncated Poisson and negative binomial distributions. 

At the beginning, we consider the abrupt cutoff only of the very high 
multiplicity tail, i.e., the case m = and iV > (n). This mimics the energy- 
momentum conservation limits. In Figs 2 and 3, it is shown that such a 
cutoff N induces oscillations of H q . The farther is the cutoff from the mean 
multiplicity, the weaker are oscillations. This quite expected result is known 



from long ago [ 12| , [15 1. It is demonstrated in Fig. 2 for (n) = 10 and different 
cutoffs at iV = 30, 40, 50. Another representation of the same result is seen 
in Fig. 3 where the constant cutoff N = 30 has been chosen for different 
(n) equal to 5, 10 and 15. The closer is the cutoff to (n), the stronger the 
low-rank moments are damped. For the faraway cutoff N = 50, the period 
of oscillations increases. This increase is larger for lower mean multiplicity 
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(see Fig. 4). At iV/(n)= const, one observes the approximate scaling of H q 
as seen in Fig. 5. 



3 Very high multiplicities 

With the advent of RHIC, LHC and TESLA we are approaching the situation 
when average multiplicities become very high and the tails of multiplicity 
distributions reach the values which are extremely large. These events with 
extremely high multiplicities at the tail of the distribution can be of a special 
interest. The tails of particular channels die out usually very fast, and a 
single channel dominates at the very tail of the distribution. Mostly soft 
particles are created in there. Thus one hopes to get the direct access to 
very low-x physics. QCD-interpretation in terms of BFKL-equation (or its 
generalization) can be attempted. Also, the hadronic densities are rather 
high in such events, and the thermodynamical approach can be applied |L6 



However, these events are rather rare and the experimental statistics is 
quite poor until now. The Poisson distribution has the tail which decreases 
mainly like an inverse factorial. According to NBD ([!]), the tail is exponen- 
tially damped with the power-increasing preexponential factor. At the same 
time, QCD predicts even somewhat slower decrease as is seen in Fig. 1 from 
the behaviour of moments. This is important for future experiments in the 
very high multiplicity region. 

To study these events within the truncated NBD with k = 2 according to 
Eqs (|i~5l), (|16D , let us choose the multiplicity interval of the constant length 
N — m = 20 and place it at various distances from the mean multiplicity 
(n) = 10 as in Fig. 2. The resulting H q are shown in Fig. 6. The most 
dramatic feature is the negative values of H 2 and the subsequent change 
of sign of H q at each q in the case when the lower cutoff m is noticeably 
larger than (n) (m/(n) > 2). This reminds of the behaviour of H q for the 
fixed multiplicity distribution and shows that the NBD-tail decreases quite 
fast so that the multiplicity m dominates in the moments of these truncated 
distributions. 

The same features are demonstrated in Figs 7 and 8 for different average 
multiplicities and different positions of the fixed window N — m = 20. In Fig. 
7, the window is rather close to (n) or even contains it inside (for (n)=15). 
Therefore H2 < only for (n)=5. In Fig. 8, it is very far from (n). Thus 
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all H 2 are negative, the more the lower is (n). Again, the sign-changing 
characteristics remind those for the fixed multiplicity distribution. 

Another possibility to study the tail of the distribution with the help of 
i/q-ratios is their variation with the varying length of the tail chosen. At 
the same mean multiplicity (ra)=10, we calculate moments for the intervals 
starting at m = 20 and ending at N = 50, 40, 30, 25. The values of H q at 
rather low ranks q = 2, 3, 4, 5 are very sensitive to the interval length as 
shown in Fig. 9. The values of H 2 vary, e.g., by the order of magnitude. 

4 Conclusions 

In connection with some experiments planned, our main concern here was 
to learn if if g -ratios can be used to judge about the behaviour of the tail of 
the multiplicity distribution. Using NBD as an example, we have shown that 
H q behave in a definite way depending on the size of the multiplicity inter- 
val chosen and on its location. Comparing the corresponding experimental 
results with NBD-predictions, one would be able to show whether the exper- 
imental distribution decrease slower (as predicted by QCD) or faster than 
NBD. 

In particular, the negative values of H 2 noted above are of special interest 
because they show directly how strong is the decrease of the tail. NBDs at 
different k values would predict different variations of H 2 with more negative 
H 2 for larger k. Also, the nature of oscillations of if g -moments at larger 
values of q reveals how steeply the tail drops down. 

Let us stress that the choice of high multiplicities for such a conclusion 
could be better than the simpleminded fit of the whole distribution. As 
one hopes, in this case there is less transitions between different channels 
of the reaction (e.g., from jets with light quarks to heavy quarks), and the 
underlying low-x dynamics can be revealed. 

Acknowledgements 

This work is supported by the RFBR grants N 00-02-16101 and 02-02- 
16779. 

Figure captions. 

Fig. I. D q — q 2 H q are shown for asymptotical QCD by circles (D® CD = 
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1) in comparison with NBD-predictions at different values of the parameter 
k=2 (diamonds), 3 (crosses), 10 (squares). The numbers near the QCD values 
show which values of k one would need to use for NBD to fit D q — 1 at the 
corresponding q. 

Fig. 2. i^-moments for NBD with (n)=10 truncated at N = 30 (dia- 
monds), 40 (crosses), 50 (squares) at m — 0. 

Fig. 3. iJg-moments for NBD with (n)=5 (diamonds), 10 (crosses), 15 
(squares) truncated at N = 30 at m = 0. 

Fig. 4. //^-moments for NBD with (n}=5 (diamonds), 10 (crosses), 15 
(squares) truncated at N = 50 at m = 0. 

Fig. 5. Approximate scaling of iJg-moments for NBD with iV/(n) =const. 

Fig. 6. ifg-moments for NBD with (n)=10 truncated at m = 10, iV = 30 
(diamonds), 20-40 (crosses), 30-50 (squares). 

Fig. 7. iJ g -moments for NBD with (n}=5 (diamonds), 10 (crosses), 15 
(squares) truncated at m = 10, iV = 30. 

Fig. 8. i/ g -moments for NBD with (n)=5 (diamonds), 10 (crosses), 15 
(squares) truncated at m = 30, iV = 50. 

Fig. 9. ^-moments for NBD with (n)=10 truncated at m = 20, iV = 50 
(diamonds), 20-40 (crosses), 20-30 (squares), 20-25 (circles). 
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